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[i] The structure of Mercury’s internal magnetic field has been determined from analysis 
of orbital Magnetometer measurements by the MESSENGER spacecraft. We identified 
the magnetic equator on 531 low-altitude and 120 high-altitude equator crossings from 
the zero in the radial cylindrical magnetic field component, B p . The low-altitude crossings 
are offset 479 ± 6 km northward, indicating an offset of the planetary dipole. The tilt of the 
magnetic pole relative to the planetary spin axis is less than 0.8°. The high-altitude 
crossings yield a northward offset of the magnetic equator of 486 ± 74 km. A field 
with only nonzero dipole and octupole coefficients also matches the low-altitude 
observations but cannot yield off-equatorial B p — 0 at radial distances greater than 3520 km. 

We compared offset dipole and other descriptions of the field with vector field observations 
below 600 km for 13 longitudinally distributed, magnetically quiet orbits. An offset 
dipole with southward directed moment of 190 nT-i?M yields root-mean-square (RMS) 
residuals below 14 nT, whereas a field with only dipole and octupole terms tuned to match 
the polar field and the low-altitude magnetic equator crossings yields RMS residuals up to 
68 nT. Attributing the residuals from the offset-dipole field to axial degree 3 and 
4 contributions we estimate that the Gauss coefficient magnitudes for the additional terms 
are less than 4% and 7%, respectively, relative to the dipole. The axial alignment and 
prominent quadrupole are consistent with a non-convecting layer above a deep dynamo in 
Mercury’s fluid outer core. 
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1. Introduction 

[2] The discovery of Mercury’s internally generated mag- 
netic field by the Mariner 10 spacecraft [Afess et al., 1974, 
1975] presented an enigma: the large-scale structure of the 
field suggested a dynamo origin, but the observed field 


'Space Department, The Johns Hopkins University Applied Physics 
Laboratory, Laurel, Maryland, USA. 

"Department of Earth and Ocean Sciences, University of British 
Columbia, Vancouver, British Columbia, Canada. 

"Planetary Science Institute, Tucson, Arizona, USA. 

4 Space Science Institute, Boulder, Colorado, USA. 
s NASA Goddard Space Flight Center, Greenbelt, Maryland, USA. 
"Department of Atmospheric, Oceanic and Space Sciences, University 
of Michigan, Ann Arbor, Michigan, USA. 

"Department of Terrestrial Magnetism, Carnegie Institution of Washington, 
Washington, D. C., USA. 

"Lamont-Doherty Earth Observatory, Columbia University, Palisades, 
New York, USA. 

"Department of Earth, Atmospheric and Planetary Sciences, Massachusetts 
Institute of Technology, Cambridge, Massachusetts, USA. 

Corresponding author: B. J. Anderson, Space Department, 

The Johns Hopkins University Applied Physics Laboratory, 1 1 100 Johns 
Hopkins Rd., MS B200-E230, Laurel, MD 20723-6099, USA. 
(brian.anderson@jhuapl.edu) 

©2012. American Geophysical Union. All Rights Reserved. 
0148-0227/12/2012JE004159 


strength was one or more orders of magnitude less than that 
predicted by scaling arguments and the fields measured at 
other planets [Stevenson, 2003], Alternative explanations 
included long-wavelength crustal remanence [e.g., Aharonson 
et al., 2004] and thermoelectric currents at the core-mantle 
boundary [e.g., Stevenson, 1987; Giampieri and Balogh, 2002]. 

[3] The recent finding that at least a layer of the core 
remains molten at present [Margot et al., 2007] favors a 
dynamo source rather than remanent magnetization for the 
origin of the field [Aharonson et al, 2004], This interpretation 
is further supported by magnetic field observations taken 
during the first two flybys of Mercury by the MErcury Sur- 
face, Space ENvironment, GEochemistry, and Ranging 
(MESSENGER) spacecraft, which revealed a dominantly 
dipolar field, with no evidence for large-scale remanent crustal 
fields [Anderson etal., 2008; Punicker et al., 2009; Uno etal., 
2009]. A variety of core structures and dynamo regimes have 
been demonstrated to yield planetary dipole moments that are 
sufficiently weak to be compatible with observations of Mer- 
cury’s magnetic field. These include thin-shell dynamos [e.g., 
Stanley et al., 2005; Takahashi and Matsushima, 2006], deep 
dynamos [Heimpel et al., 2005], deep dynamos enclosed by a 
stably stratified, electrically conductive layer [Christensen, 
2006; Wicht et al, 2007; Christensen and Wicht, 2008; 
Manglik et al., 2010], induction feedback on the converting 
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portions of the core from magnetopause currents [Glassmeier 
et al, 2007a, 2007b; Gomez-Perez and Solomon, 2010; 
Gomez-Perez and Wicht, 2010; Heyner et al., 201 la, 201 lb], 
and precipitation of solid iron in radial zones within a liquid 
outer core [ Vilim et al., 2010]. 

[ 4 ] Recent orbital observations from MESSENGER con- 
firm that the planetary field can be represented as a nearly 
axially aligned dipole. Elowever, not only is the dipole 
moment weak, but it is also offset by 0.19 R M (where R M is 
Mercury’s mean radius, 2440 km) to the north. This offset is 
equivalent to axial dipole, gi 0 , and quadrupole, g 2 0 , terms in a 
spherical harmonic representation in the ratio g 2 Jg\o = 0.4 
[Anderson et al., 2011], Such a result is consistent with prior 
modeling from flyby observations [Uno et al., 2009; Alexeev 
et al., 2010] and helps clarify why the MESSENGER equa- 
torial flybys, inbound from the tail dusk side, all passed 
through the southern magnetotail lobe [Anderson et al, 2010], 

[5] The combination of a large quadrupole term with close 
alignment of the weak field to the planetary rotation axis 
presents a challenge to existing dynamo models. Downward 
continuation of an axial field with the inferred gnjgio ratio 
yields a field at the core-mantle boundary that is still dipole- 
dominated but with a large quadrupole term. Conventional 
dynamos, in which the top of the convecting layer coincides 
with the top of the core, have difficulty in producing such a 
field. Multipolar dynamos can produce a weak field [e.g., 
Olson and Christensen, 2006], but such fields are not dipole 
dominated. Thin-shell dynamos can produce weak fields but 
show a dipole offset much less than observed at Mercury 
[Stanley et al., 2005], or equivalently no clear dominance of 
harmonic degree 2 (i.e., quadrupole) over higher degrees in 
the power spectra [Heimpel et al., 2005; Takahashi and 
Matsushima, 2006], The feedback-dynamo process yields a 
planetary field that is axisymmetric and dominated by odd 
hannonics [Heyner et al., 201 lb], and the influence of iron 
snow layers also produces a field with a large octupole 
(degree 3) contribution and comparatively weak quadrupole 
structure [Vilim et al., 2010], In another class of models, sub- 
adiabatic heat flow at the core-mantle boundary produces a 
stably stratified layer. This stable, electrically conductive 
layer filters out the spherical harmonic terms of high degree 
and order that are strong in the field at the top of the dynamo 
region, but that also vary the most in time. In these models the 
field at the planetary surface is dominated by the axial dipole 
and quadrupole terms that are the most slowly time-varying 
components [Stevenson, 1982; Christensen, 2006; Wicht 
et al., 2007; Christensen and Wicht, 2008]. In addition, 
zonal flows set up in the stably stratified layer preferentially 
reduce non-axisymmetric, low-degree terms [Christensen 
and Wicht, 2008], yielding a surface field that is weak, is 
dipole dominated, has a substantial quadrupole component, 
and is axisymmetric. Although these models are promising, 
the stable layer thickness, its dynamics (e.g., the role of 
double diffusion), and the resulting heat transport have been 
shown to affect markedly the observable field structure 
[Stanley and Mohammadi, 2008; Manglik et al, 2010]. 

[6] The implications of the offset dipole result therefore 
warrant a critical examination of the initial findings from 
Mercury orbital data. The Anderson et al. [2011] analysis 
identified the magnetic equator from the spacecraft position 
at the point where the radial component of the magnetic field 
(B p ) vanished in a cylindrical coordinate system in which the 


cylinder axis, or z-axis, is aligned with the planetary spin 
axis, i.e., where the magnetic field was parallel to the spin 
axis. Interpretation of these observations is not unambigu- 
ous, however, because it is possible for more complex field 
geometries not indicative of an offset dipole to yield surfaces 
on which B p = 0. For example, a field with a dipole and large 
octupole terms but a small quadrupole term, as predicted by 
the induction feedback mechanism of Heyner et al. [201 la, 
2011b], will yield three surfaces on which B p = 0, and it is 
possible that only the most northern of these might be 
detectable with a spacecraft in an orbit as eccentric as that of 
MESSENGER. 

[7] In this paper we examine the MESSENGER magnetic 
field data from Mercury orbit in detail to assess the 
uniqueness of the B p = 0 signature as an indicator of the 
magnetic equator. We also present initial analyses of higher- 
degree structure in the axial field to more firmly distinguish 
among candidate models for the planetary field generation. 

[8] In Section 2, we present an in-depth analysis of the 
magnetic equator signatures, including distant crossings of 
the magnetotail current sheet. The possibility that the 
observed signatures reflect centered, axially aligned, odd 
harmonic fields is assessed quantitatively. In Section 3, the 
analysis is expanded to a comparison of the observed vector 
fields over all longitudes against the most likely candidates 
for the planetary field to quantify the degree of correspon- 
dence between the global vector magnetic field and different 
hypothetical global fields. The results are summarized and 
discussed in Section 4. 

2. Magnetic Equator Analysis 

[9] Here we present a more extended analysis of the B p = 0 
signature than in our earlier report [Anderson et al, 201 1] by 
considering both the descending-node, lower-altitude equator 
crossings as well as the more distant ascending-node equator 
crossings. The magnetic equator determination is particularly 
valuable for analysis of MESSENGER orbital observations. 
Using only multivariable, nonlinear inversions to solve for 
Mercury’s internal field is problematic for several reasons 
[cf. Johnson et al., 2012], First, high plasma pressures and 
corresponding magnetization currents are prevalent below 
0.5 R m altitude and within ~45° latitude of the equator and 
cause depressions in the magnetic field strength from tens of 
percent to over 90% [Anderson et al, 2010; Korth et al., 
2011, 2012; Winslow et al., 2012]. Second, external fields 
make contributions comparable to the internal field even at 
the planetary surface, so their effects must be quantitatively 
included [e.g., Alexeev et al., 2008, 2010]. Because plasma 
pressure locations and densities vary from orbit to orbit [e.g., 
Korth et al., 2012], however, the errors in the external field 
corrections for individual orbits can be large. Moreover, one 
cannot solve for the external field using the potential for- 
malism because the sampled region is not free of currents. 
Finally, the MESSENGER orbit is eccentric and there is a 
northern hemispheric bias in low-altitude sampling, leading 
to high covariance among terms of low degree and order in 
the inversions [Korth et al., 2004]. Restricting data to lati- 
tudes poleward of 45°N exacerbates the covariance problem. 

[ 10 ] We therefore opted to determine the magnetic equator 
first and use this as a constraint on the inversions. For 
Mercury, the flyby observations showed that the dipole field 
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Figure 1. MESSENGER magnetic field and trajectory data for the Mercury magnetospheric transit on 
28 March 2011, day of year (DOY) 087, from 01:30 to 03:30 UTC. (top) The 1-s averaged magnetic field 
data in MSO cylindrical coordinates (B p , B B z ) together with the field magnitude, B T , and its negative, 
— B t . (bottom) Spacecraft position in MSO coordinates: latitude, longitude, and altitude above the plane- 
tary surface. Longitude is measured positive eastward relative to the planet-Sun line such that —90° is 
dawn, 0° is noon, 90° is dusk, and 180° (—180°) is midnight. The bow shock and magnetopause bound- 
aries are indicated by vertical dashed lines labeled SK and MP, respectively. The descending-node plan- 
etary equator crossing, 0° geographic latitude, is also indicated. 


is close to axially aligned [Afess et al., 1975; Anderson et al., 
2008], and this geometry yields a particularly well-defined 
magnetic equator. If the planetary moment is axially aligned 
and perpendicular to the solar wind flow direction, the 
external currents (magnetopause and cross-tail currents) are 
centered on the magnetic equator and do not influence its 
location. Thus, for Mercury, the magnetic equator can be 
determined without making any corrections for the external 
fields and is an accurate independent constraint for both the 
external field description and inversion analyses [cf. Johnson 
et al., 2012], 

[l l] For our analysis, we use Mercury solar orbital (MSO) 
coordinates, for which the Cartesian (X, Y, Z) components 
are as follows: Z is positive along the planetary spin axis, X 
is positive toward the Sun, and Y completes the right-hand 
system, positive duskward. Data are also analyzed in MSO 
cylindrical coordinates ( p , </>, Z), where the unit vectors, n, 
are defined as follows: n^, is positive eastward in the plane 
parallel to the planetary equator, n z is northward, and n p 
completes the right-hand system such that n p = n^, x n z . 

[ 12 ] The MESSENGER spacecraft [ Solomon et al., 2007] 
entered orbit about Mercury on 18 March 2011. The initial 
spacecraft orbit had a periapsis altitude of 200 km, 82.5° 
inclination, 15,300 km apoapsis altitude, and 12-h period 
[cf. Zurbuchen et al., 2011], Data in this study are from the 
MESSENGER Magnetometer [Anderson et al., 2007] 
from 23 March 2011, shortly after orbit insertion, through 
2 March 2012. For this date range the ascending and des- 
cending orbit segments crossed the planetary equator at 
altitudes from 1000 to 1500 km, and from 3500 to 5000 km, 
respectively. 


2.1. Descending-Node Equator Crossings 

[13] Magnetic field and trajectory data are shown in 
Figure 1 for the first magnetospheric transit on 28 March 
2011. The bow shock and magnetopause crossings identified 
from the rotations in the magnetic field are indicated along 
with the planetary equator crossing. The spacecraft was 
inbound in the afternoon near 45° MSO longitude, crossed the 
descending orbit node in the early morning near — 135° MSO 
longitude, and exited the magnetosphere 18 min later. The 
maximum field magnitude occurred slightly before the mini- 
mum altitude, consistent with the contributions from the 
magnetopause and tail current systems. 

[14] On the descending orbit segment near 02:30 UTC the 
B p component reversed sign from negative to positive sev- 
eral minutes prior to and north of the planetary equator 
crossing. When B p = 0 nT near 02:27:30 UTC, the field was 
nearly purely northward with B^ ~ 0 nT. An expanded-time 
view of the equator crossing in Figure 2 shows the same 
magnetic field and spacecraft latitude data but instead of 
planetocentric distance gives the angle between n p and the 
magnetic field defined as cos -1 (Bp/Bj), where Bj is the total 
magnetic field magnitude. This angle changed from near 
160°, reflecting an inward-directed field, passed through 90° 
at the magnetic equator, and then decreased below 30° as the 
spacecraft passed into the southern magnetotail lobe. The 
magnetic equator crossing preceded the planetary equator by 
~ 150 s, and at the planetary equator B p was ~30 nT. 

[15] To quantify the spacecraft position when B p passes 
through zero, we consider Z MSO as a function of B p as shown 
in Figure 3 for four sequential descending-node equator 
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Figure 2. Expanded plot of MESSENGER Magnetometer data from Figure 1, 28 March 2011, focusing 
on the planetary equator crossing near 02:30 UTC. (top) The 1-s averaged magnetic field in MSO cylin- 
drical coordinates in the same format as the top panel in Figure 1. (bottom) Angle between the magnetic 
field and n p , cos ~ l (B p /B T ), and spacecraft MSO latitude. The vertical dashed lines indicate times when the 
spacecraft passed through the magnetic equator (close-spaced dashes) and the planetary equatorial plane 
(wide-spaced dashes). The B p = 0 point occurs at a planetary latitude of ~10° and ~ 1 50 s prior to the 
equator crossing. 


2011:086 02:18:50 , 2011:086 14:22:52 




2011:087 02:27:23 



2011:087 14:32:01 



Figure 3. Spacecraft Z MSO position plotted versus B p for four segments of data near the descending-node 
equator crossings for four consecutive orbits from 27 and 28 March 201 1 (DOY 86 and 87). Observations 
are plotted as plus symbols connected by solid thin black lines, and quadratic fits are shown by the thick 
gray lines. The Z p0 values for each pass obtained by the fits are indicated in each panel; uncertainties in the 
fit intercepts are three-standard-error estimates. 
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Figure 4. Magnetic equator crossing points, Z p 0 , versus (a) MSO local time and (b) planetary longitude 
derived from all MESSENGER Magnetometer data from Mercury orbit between 24 March 2011 and 2 
March 2012. Vertical error bars show three-standard-error estimates in Z p o- The lack of values in the morn- 
ing hours results from the fact that MESSENGER encounters the magnetopause before crossing the mag- 
netic equator. Note that the average apparent local time of the solar wind ram direction is near 1 1 h MSO, 
so that the times of sparse data are nearly centered on the solar wind ram direction. The systematic vari- 
ation in Z p0 with either local time or longitude has an amplitude of ~50 km (100 km peak to peak), about 
10% of the offset. 


crossings on 27 and 28 March 2011. The example from 
Figures 1 and 2 is in Figure 3c. These cases also show 
quadratic fits to Z mso . The B p = 0 intercept of each fit gives 
an estimate of Z at the magnetic equator crossing for that 
pass, denoted as Z p0 . The standard error estimate of the 
intercept is used as a measure of confidence in the Z p0 . Note 
that near the B p = 0 point we do not observe dB p /dZ ~ 0. 
This is an important feature to be compared with the forward 
modeling of dipole-octupole models considered in Section 4. 

[16] Corresponding fits were performed to all of the 
descending-node equator crossing data from 24 March 2011 
to 2 March 2012. The results for the 531 near-planet equator 
crossings identified are shown in Figure 4. When the des- 
cending node was on the dayside, the spacecraft often 
encountered the magnetopause either before or nearly coinci- 
dent with the B p = 0 point, so Z p o values were not obtained in 
these cases. The Z p0 data suggest a possible ~50 km amplitude 
variation such that Z p0 is greater at noon than at midnight. The 
average radial spacecraft distance corresponding to these 
crossings was ~3400 km, so the variation in Z p0 corresponds 
to a peak-to-peak latitude variation of ~1.5°. The mean Z p o, 
the mean standard error, and the average of Z p0 weighted by 
the inverse of the standard errors are given in Table 1 together 
with the altitude range of the points. The mean and weighted 
mean Z p0 values are not significantly different, but we consider 
the weighted mean to be slightly more reliable and adopt the 
value 479 ± 6 km as the offset value from these data. 

[ 17 ] With this value for the offset we next assess the tilt in 
the field relative to the planetary rotation axis. To quantify 
the tilt of the offset dipole we transform the position and 


magnetic field data into a body-fixed coordinate system 
offset 479 km to the north along the planetary rotation axis, 
which we term the Mercury MAGnetic (MMAG) coordinate 
system, and we evaluate the latitude for which the radial 
field component, R r _ M MAG> in this system vanishes. The 
latitudes of the 5,-mmag = 0 points give the magnetic dip 
latitude of each crossing in the MMAG system. The results 
of this analysis are shown in Figure 5. To determine the 
plane that contains the dip equator, we fit these data with a 
sine function, and the amplitude of this fit yields a tilt of 

Table 1 . Average Magnetic Offset Zpo Values, Uncertainties, and p 
Range and Mean Values for Descending (Near) and Ascending (Far) 
Equator Crossings Derived From MESSENGER Magnetometer Data 
From 24 March 2011 Through 2 March 2012 a 



Descending (Near) 

Ascending (Far) 

Units 

Number of cases 

531 

120 

n/a 

Mean Z f)0 =t 3 a 

All ± 7 

420 ± 86 

km 

Mean 3 a 

24 

87 

km 

Weighted mean Z p o ± 3a 

479 ± 6 

486 ± 74 

km 

Mean p MSO ± sdev 

3410 ± 110 

6110 ± 470 

km 

Pminj Pmax 

3150, 3720 

5240, 7560 

km 


a The uncertainties in the mean Z p0 values are three standard errors of 
the mean (3cr)- The mean 3cr values are the averages of the 3<r standard errors 
in the Z p o fit value for each equator crossing. The weighted mean Z /)0 is the 
average Z p 0 weighted by the inverse of the fit value standard errors. For the 
Pmso statistics, sdev is the population standard deviation, whereas p m j n and 
Pmax ar e the minimum and maximum p MSO values in each population. 
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Figure 5. Dip equator latitude plotted versus planetary longitude. Data are shown by X symbols and 
error bars are three standard errors. The gray curve shows a sine fit constrained to have a period of 
360° and corresponds to the best fit plane that represents the tilt of the magnetic equator relative to the 
planetary rotation axis. The resulting tilt is 0.8° with an azimuth (for the southern hemisphere pole) of 
—77° longitude. Because the data do not exhibit a strong 360° period, this tilt is considered an upper limit. 


0.8°. Since the tilt signal is weak relative to the scatter in the 
data, we consider this result to be an upper limit, and we do 
not regard the longitude of the tilt as sufficiently constrained 
to report a value for it. 

2.2. Ascending-Node Equator Crossings 

[is] Here we present an analysis of the more distant 
crossings, which occurred from 3500 to 5000 km altitude in 
the magnetotail. If the near-planet equator crossings corre- 
spond to a taie offset of the dipole, we expect that the more 
distant crossings would display a similar offset relative to the 
planetary equator. However, if the lower-altitude crossings 
were due to complex higher-degree structure in the field, 


these higher-degree terms should fall off more rapidly than 
the dipole and the higher-altitude magnetic equator would be 
shifted toward the planetary equator. 

[19] An example ascending-node equator crossing on 
24 November 2011 from 08:00 to 11:30 UTC is shown in 
Figure 6. The distant equator crossings occur in the mag- 
netotail where the spacecraft transitions from the southern 
magnetotail lobe, within which the magnetic field is nomi- 
nally anti-sunward and B ~ _3 T n p , to the northern magnetotail 
lobe where the field is nominally sunward and B ~ — B T n p . 
Coincident with the reversal in B p is a depression in the total 
field indicative of the plasma pressure in the current sheet 
[Korth et al., 2011, 2012]. In this particular instance the 


MESSENGER Magnetometer 



Figure 6. MESSENGER magnetic field data for the Mercury magnetospheric transit on 24 November 2011 
(DOY 328) from 08:00 to 1 1 :30 UTC in the same format as Figure 1 , with the exception that the inbound bow 
shock crossing is omitted for clarity. The inbound bow shock crossing occurred at 05:45 UTC. The ascend- 
ing-node planetary equator crossing, 0° latitude, is indicated near 09:23 UTC. The spacecraft was inbound in 
the late evening and crossed the ascending orbit node near midnight, at ~ 1 70° MSO longitude. The outbound 
magnetopause was encountered prior to crossing the descending node equator. 
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Figure 7. Expanded plot of MESSENGER Magnetometer data from Figure 6, 24 November 2011, 
focusing on the planetary equator crossing near 09:23 UTC in the same format as Figure 2. 


current-sheet crossing started several minutes after the space- 
craft crossed the magnetic equator, and the center of the cur- 
rent sheet was crossed ~5 min after the planetary equator. The 
maximum magnetic field occurred slightly after the minimum 
altitude, consistent with the dayside enhancement of the polar 
field due to magnetopause currents. For the example in 
Figure 6, the outbound magnetopause was encountered prior 
to the descending-node equator crossing. The outbound shock 
crossing occurred at 10:40 UTC, although variability, which we 
attribute to foreshock dynamics, was evident until 1 1 :00 UTC. 

[20] An expanded view of the ascending-node equator and 
current sheet crossings on this orbit is shown in Figure 7. 
The depressions in B T coincide with rotations in the field 
from the southern tail lobe (B ~ B p n p ) to the northern lobe 
(B rxj — B p n p ) as the spacecraft passed through the magneto- 
tail plasma sheet and cross-tail current sheet. The current 
sheet crossing occurred over a span of ~7 min starting 
slightly more than 3 min after the planetary equator crossing. 
The azimuthal component in the field remained small 
throughout the crossing. In contrast to the lower-altitude 
equator passes, the passes through the deeper tail current 
sheet often exhibited multiple reversals in B p , which we 
attribute to dynamics in the tail related to north-south motion 
of the tail, magnetic reconnection dynamics, and turbulence 
in the current sheet. 

[21] The higher- altitude crossings were analyzed in the 
same manner as the lower-altitude crossings to determine a 
best fit Z p0 for each crossing. Plots of Z MSO versus B p are 
shown in Figure 8 for four sequential distant equator crossings 
on 24 and 25 November 2011. The greater variability in B p for 
the distant tail crossings is evident, reflecting repeated passage 
of the current sheet back and forth over the spacecraft. The fit 
intercepts yield best estimates for the average equator for each 
crossing. The standard errors of the intercepts are correspond- 
ingly higher for these distant tail crossings. 

[22] To compare the descending (near) and ascending (far) 
equator crossings we plot Z^o versus Pmso for all of the cross- 
ings in Figure 9. Due to the evolution of the MESSENGER 


spacecraft orbit, successive deep tail (or “long-eclipse”) seasons 
occurred for trajectories progressively closer to the planet. 
There are only six equator crossings determined during the first, 
most distant, deep tail season, because the Magnetometer was 
powered off for most of the first of these seasons. The deep tail 
seasons correspond to times of longest passage of the spacecraft 
through the planetary shadow, times during which management 
of spacecraft power was needed. For the first tail season, the 
Magnetometer instrument was powered off to safeguard the 
spacecraft while the behavior of the power system in Mercury 
orbit during long-eclipse seasons was characterized. Fortu- 
nately, the power margins allowed Magnetometer operation 
through successive long-eclipse seasons. 

[23] The scatter in the more distant tail crossings is consid- 
erable and is larger than the three-standard-error uncertainties, 
so the variability must be real. Nonetheless, the crossings are 
not centered on the planetary equator. The statistics for the 
distant tail crossings are given in Table 1. The mean and 
weighted mean values for the near and far equator crossings 
are in agreement, indicating that the cross-tail current above 
3000 km altitude is offset to the north. This result is consistent 
with expectations for an offset dipolar planetary field that 
organizes the global-scale external current systems including 
the tail magnetopause and cross-tail currents such that the 
mean current sheet position corresponds to the dipole mag- 
netic equator [cf. Alexeev et ai, 2008, 2010]. 

[24] The scatter in Z p o for the distant equator crossings is 
larger than that for the near tail crossings, and we briefly 
consider the source of this variation. We first examined the 
IMF fix-Mso component derived from MESSENGER obser- 
vations before and after each magnetospheric transit. As 
Winslow et al. [2012] have shown, the sign of the IMF Rx-mso 
component was generally consistent on the inbound and out- 
bound MESSENGER spacecraft orbits, so this component 
should reliably indicate the IMF sector imposed on Mercury’s 
magnetosphere at the time of the equator crossings. We found 
no correlation between Z p o and IMF Rx-mso, indicating that 
this quantity is not governing the scatter in Figure 9. 
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Figure 8. Spacecraft Z MSO position plotted versus B p for four segments of data near the ascending equa- 
tor crossings from four consecutive Mercury orbits on 24 and 25 November 201 1 (DOY 328 and 329) in 
the same format as Figure 3. 


[25] We then divided the distant equator crossing events into 
three ranges of Z p0 with nearly equal numbers of events in each 
range, and we evaluated the average magnetic fields and posi- 
tions as functions of time relative to the time of the current sheet 
crossing. The superposed epoch averages as functions of time 
relative to the fit B p = 0 point are shown in Figure 10. The 
average \B (p \ was very small, so we show only B p and B z 
together with the total field, B T . The angle relative to northward 
of the field projected onto the p-Z plane, 9 P . Z = tan ^ l (Bp/B z ), is 
also shown together with the difference between 6 p _ z for the 
southernmost and northernmost events, 89 p _ z = 9 p _z _ s — 9 p _ Zt N . 
The number of events, the average, minimum and maximum 
Zpo, and the mean 9 p _ z for the low, medium, and high Z p o 
ranges used are given in Table 2. 

[26] The Z component in the southern lobe, before the 
current sheet crossing, is consistently more positive for the 
northernmost than for the southernmost events. For the most 
southern Z p0 events, the magnetic field is tilted by 5° to 10° 
toward more positive 9 p _ z than for the most northern Z p0 
events. The mean 89 p . z values for the southern lobe from 
1500 to 1000 s and 1000 to 500 s before the current sheet 
crossing are 6.4° ± 0.2° and 5.3° ± 0.2°, and that for the 
northern lobe from 500 to 1000 s after the current sheet 
crossing is 6.06° ± 0.12°. These values are consistent with 
an overall southward tilt of the magnetotail for the southern 
Z p o events relative to the northern Z p o events and indicate 
that the spread in the Z^o values is due to tilting of the tail 
northward or southward. 


[27] To assess the source of the inferred north-south tail 
tilt, we consider the north-south excursions in solar wind 
velocity. Although MESSENGER does not measure the solar 
wind velocity, data from the Helios 1 and 2 spacecraft pro- 
vide information regarding the excursions in solar wind 



Figure 9. Magnetic equator crossing points, Z p0 , plotted 
verses Pmso for all crossings identified in data from 24 March 
2011 through 2 March 2012. Vertical error bars indicate three 
standard errors for each Z p0 crossing point. Black curves 
show average trajectories for the four deep-tail seasons, and 
the thick light-gray curves show trajectory segments for the 
innermost and outermost orbits from which the lower-altitude 
equator crossings were identified. 


8 of 17 


E00L12 


ANDERSON ET AL.: MERCURY’S PLANETARY MAGNETIC FIELD 


E00L12 



Figure 10. Superposed epoch analysis of magnetic field geometry for distant tail crossings. From top to 
bottom: total field |B|, B p , B z , and — |B| for the northernmost (N) distant tail crossings; same quantities for 
the southernmost (S) distant tail crossings; 9 p _ z , the tilt of the field in the p-Z plane from northward, with 
positive angles for positive p\ the difference between 6 p _ z for southernmost and northernmost distant tail 
crossings. 


velocity relative to the ecliptic plane at Mercury orbital dis- 
tances. Plasma measurements [Rosenbauer et al., 1977] in the 
slow solar wind from the Helios 1 and Helios 2 spacecraft were 
used to determine the angular deviation of the solar wind flow 
velocity from the radial direction [cf. Borovsky, 2012], Using 
40 s averages of the plasma flow velocity, the north-south 
deviation of flow was measured as ^north-south = sin _1 (fsw-z/ 
fswX where V S w is the solar wind speed and V S w-z is the 
component of the solar wind velocity perpendicular to the radial 
direction and out of the ecliptic, positive northward. In 
Figure 11, the distributions of 0 n oith-south are shown for two 
heliocentric distance ranges: 0.9-0.98 AU (blue curve) and 
0.29-0.35 AU (red curve). The mean value and root-mean- 
square (RMS) deviation of the angular deviations are denoted in 
the figure. In the distance range 0.29-0.35 AU from the Sun, 
the measured deviations have an RMS spread of 3.9°, whereas 
near 1 AU the distribution is slightly less broad with an RMS of 


Table 2. Statistics for Three Z p0 Ranges of Distant Tail Magnetic 
Equator Crossings 3 


Z p0 Range 

Low 

Medium 

High 

Units 

Number of cases 

40 

42 

40 

n/a 

Z p o avg (min, max) 

72 

493 

792 

km 


(-1054, 338) 

(356, 604) 

(614, 1320) 


6 P . z southern lobe 

82.3 ± 0.3 

80.3 ± 0.2 

77.0 ± 0.2 

deg 

Op . z northern lobe 

-71.2 ± 0.2 

— 77.1± 0.2 

-77.3 ± 0.2 

deg 


“From top to bottom the rows give: the number of cases in each range; the 
average (avg), minimum (min), and maximum (max) Z p0 ; and the average 
8 p _z = tan - ( B f J By) for the southern and northern tail lobe intervals. The 
9 p _ z averages were taken 500 to 1000 s before the magnetic equator crossing 
for the southern lobe and 500 to 1 000 s after the magnetic equator crossing 
for the northern lobe. Uncertainties for 9 p _ z are three standard errors of the 
means. 
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Figure 11. Angle of solar wind velocity relative to the 
ecliptic plane, # no rth-southi derived from Helios 1 and 2 plasma 
observations, where positive angles are northward. Red curve 
shows the distribution for heliocentric distances from 0.29 to 
0.35 AU, relevant to Mercury, and the blue curve shows the 
distribution from 0.9 to 0.98 AU heliocentric distance, com- 
parable to Earth’s orbit. The mean and root-mean-square 
^north-south values are indicated in the figure. 


3.0°. The expected spread in north-south solar wind velocity is 
therefore consistent with the difference in magnetotail field 
observed for the southern and northern MESSENGER tail 
crossings. The magnetotail orientation will follow the solar 
wind flow direction, so we attribute the tail direction variability 
and the north-south variations in Z p0 for the distant tail cross- 
ings to north-south excursions of the solar wind velocity. 

2.3. Quantitative Assessment of Dipole-Octupole 
Field Geometry 

[ 28 ] To quantify the implication of the consistent Z M so 
offset of the magnetic equator signatures in the near and distant 
tail, we consider the structure of more complex axially aligned, 
centered magnetic fields. An axially aligned and centered 
dipole-octupole field, that is, with terms for degrees 1 and 3 in 
a spherical harmonic expansion, yields three surfaces on 
which B p = 0, one of which is the planetary equator and the 
other two are symmetrically displaced to the north and south of 
the planetary equator. The ratio of the octupole to the dipole 
term falls off with distance r from the planet center of mass as 
1/r 2 , so these surfaces should converge to the planetary equa- 
tor with distance, and beyond some distance only the dipole 
B p = 0 surface should remain. The ratio of average Pmso 
values for the near and distant equator crossing populations is 
1.8, so any octupole contribution is a factor of 3.2 less at the 
location of the ascending-node equator crossings than at the 
descending-node equator crossings. Also, north-south sym- 
metric tenns of odd higher spherical harmonic degree, 5 or 
greater, decrease at least as fast as Hr relative to the octupole 


term, so the structure at MESSENGER orbit distances should 
be dominated by the dipole-octupole structure. 

[ 29 ] We therefore consider a field with only gi 0 and g 30 
terms in the spherical harmonic expansion of the magnetic 
potential, 1 // 

00 \ J+i 1 

V = r mY^ ( — ) ( Slm C0S ( m( P) + h im sin (m<p))P"'( cos 9), 

(1) 

where (r, 8, cp) is the planetocentric position in spherical 
coordinates and the P ” are the Schmidt semi-normalized 
Legendre functions. Including only the gi 0 and g 30 terms and 
writing the associated Legendre functions explicitly one finds 

5 p = fy^w| 3 ^o + ^-( 35c °s 2 6'- 1 5)g3o| (2) 

where r M = r/R M . Setting G 30 = g 3 ofyio it immediately fol- 
lows that the off-equator zeroes in B p occur for 

(3) 

and that multiple zeroes in B p occur only for radial distances 
less than a critical distance r c , i.e., where 


f'M < pcA'm 



(4) 


To yield B p = 0 at the location of the lower-altitude equator 
crossing, p M so = 3409 km and Z MSO = 477 km, requires 
G 30 = 0.546 as obtained from (3). For this ratio, the maximum 
p at which an off-equatorial zero in B p occurs is 3522 km, 
well below the minimum altitude of the distant equator 
crossings reported here. 

[30] The curve for which B p = 0 north of the planetary 
equator in a dipole-octupole field consistent with the average 
lower-altitude B p = 0 point is shown in Figure 12 together 
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Figure 12. Plot of the average magnetic equator for the 
near-tail and deep-tail seasons (markers with error bars) 
together with Z MSO points where a dipole-octupole field 
gives off-equatorial B p = 0 (solid gray curve). Vertical error 
bars are 3 standard errors of the mean for each population of 
Z P Q values. Horizontal error bars show the Pmso range of the 
observations. The dipole-octupole field was chosen to yield 
B p = 0 at the mean location of the near-tail observations. 
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Table 3. Relationships Between Gauss Coefficients in Planet- 
Centered (Unprimed Quantities) and Axially Displaced (Primed 
Quantities) Coordinates in Terms of the Axial Offset, A z 0 , and 
Offset Magnetic Dipole M = ( M , 0 m , (f> m ) 


Degree 

Planetocentric Coefficients 

Offset Coefficients 

1 

© — — 
"bo "bo ^ 
II II II 
© — — 
bo bo 

g' 10 = Mcos(0 m ) 
g'u = Msin(0 m )cos(<j> m ) 
h'u = Msin(6» m )sin(<> m ) 

2 

S 20 ~ 2Azo g' 10 
gn =g'nAz 0 A/3 
h 21 = h' 11 AZ 0 /V 3 
g22 = ^22 = 0 

n/a 

3 

£30 = 3Azjjg'i 0 

g31 = g'llAro/ 3/6 

h 3l = /iji Acj/vVi 

g3m = h 3 m = 0 for »1 >1 

n/a 

4 

§40 = 4Arj5g' 10 

n/a 


with the mean equator determinations from the low- and 
high-altitude populations of equator crossings. The horizon- 
tal lines for the near and distant magnetic equator crossings 
show the range of Pmso f° r the observations. It is clear that 
the distant B p = 0 observations cannot be explained under a 
dipole-octupole hypothesis, and the distant tail equator 
crossings are strong evidence for the offset in the planetary 
dipole field. 

2.4. Offset Dipole and Gauss Coefficient Descriptions 

[31] Representing the internal planetary field as an offset 
dipole is equivalent to a description in Mercury body fixed 
coordinates that is a series expansion in spherical harmonics 
with the contribution of each spherical hannonic tenn 
described by the Gauss coefficient. Table 3 gives the rela- 
tionship between Gauss coefficients in the body-fixed frame 
and a dipole with moment M directed with a polar angle 9 m 
and azimuth cf> m , offset along the Z axis by Az 0 . The Gauss 
coefficients of the offset dipole in a coordinate system cen- 
tered on the dipole, denoted by primes, are just the degree- 1 
terms, g' 10 , g' u , and h' u . Bartels [1936] discussed the 
relationship between Gauss coefficients evaluated in two 
coordinate systems for which the axes are aligned but the 
origins are displaced. Here we consider only an axial dis- 
placement, Az 0 , along the z axis, expressed as a fraction of 
the mean planetary radius, and Table 3 gives the Gauss 
coefficients in the planet-centered system in terms of, Az 0 , 
g' 10, g' 11, and h' n . The offset dipole yields Gauss coeffi- 
cients of all degrees in the planet-centered system. The 
degree-one (dipole) coefficients in the offset and planet- 
centered systems are identical. Degree / terms are propor- 
tional to the product of the dipole term and A z' 0 . Only the 
m = 1 non-axial tenns are nonzero, because a tilted dipole 
has only order-one azimuthal structure. 

[ 32 ] The results for the Mercury offset dipole obtained above 
are summarized in Table 4, which gives the offset dipole 
moment and orientation together with the corresponding values 
for the Gauss coefficients in the planet-centered system. 
Uncertainties are three-standard-error estimates. Since we 
derived only an upper limit on the dipole tilt of 0.8°, 
corresponding to 9 m > 179.2°, we do not report a value for <f> m . 


The tilt upper limit corresponds to upper limits on the magni- 
tude of the rn = 1 off-axis terms. We report the off-axis upper 
limits only through 1 = 3 because they are already negligible for 
1 = 2. Importantly, the large offset of the dipole along the 
rotation axis has substantial axial (m = 0) terms in the body- 
fixed spherical hannonic expansion up to and including the 
degree-4 coefficient, g 40 . We emphasize that the values and 
limits in Table 4 are those implied by the offset dipole field. 
By superposition, contributions to coefficients of higher degree 
and order derived by subtracting from the observed planetary 
field the field of the offset dipole, may be summed with the 
values given here. We denote such residual contributions with 
an asterisk (*). 

3. Forward-Model Comparison of Alternative 
Internal Fields 

[ 33 ] The analysis above was based solely on the equator 
crossing points identified from B p = 0. Here we consider data 
over as much of the orbit as possible to measure the agreement 
or disagreement between different models for the internal 
field. We use forward modeling to identify the dominant 
characteristics of the internal field. Because the external and 


Table 4. Properties of Mercury’s Internal Planetary Magnetic 
Field 3 


Quantity 

Value 

A zq 

479 ± 6 km 

AzJR m 

0.196 ± 0.002 

M 

190 ± 10 nT-i?M 

e m 

>179.2° 

0m 

n/a b 

gio 

-190 ± 10 nT 

Vfen + AiiVgio 

<0.014 

Vfe? 1 + An) 

<2.7 nT 

gialgxo 

0.392 ± 0.010 

§20 

-74.6 ± 4.0 nT 

^(§21 + hh)/g w 

<0.0016 

^(g2 1 + ^2l) 

<0.3 nT 

g22, ^22 

0 C 

§3()lg 1 0 

0.116 ± 0.009 

g30 

-22.0 ± 1.3 nT 

Vfeli + hh)/g w 

<0.0016 

V(g3i + hi i) 

<0.3 nT 

§ 3 m. h im for m >1 

0° 

Ig30*/gl0l 

<0.04 

§40/g\0 

0.030 ± 0.005 

g40 

-5.7 ± 0.2 nT 

lg40*/glol 

<0.07 


a Values for the dipole moment, axial offset, Ar 0 , and magnetic dipole 
M = (M, 0 m , o m ) as well as the corresponding Gauss coefficients are 
given in planet-centered coordinates. Uncertainties are 3cr standard errors. 
The quantities g 30 * and g 40 * are the contributions to the axial degree 3 
and 4 coefficients in addition to those implied by the offset dipole. 

b Reported coefficients are for zero tilt, 0 m = 180°, for which the tilt 
longitude is not defined. The upper limit on the tilt of 0.8° was not judged 
sufficiently robust to warrant evaluation of (f> m . 

c An axially offset dipole does not yield nonzero tenns for orders greater 
than 1 . Signals in residual fields from the axial offset dipole corresponding 
to non-axial coefficients of order higher than 1 are not estimated here, since 
the longitudinal structure in residuals appears to be ordered better in 
magnetospheric coordinates than in body-fixed coordinates. 
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Table 5. Representative Magnetic Equator Crossings With Low 
Standard Errors in •Zpo for Which the Magnetospheric Transits 
Were Used for Forward Model Comparisons 3 


Date (DOY) 

B p = 0 Time 
(hhmm) 

Zpo 

(km) 

3 a Zpo 
(km) 

MSO 
LT (h) 

Longitude 

(°E) 

26 Mar 2011 (85) 

1414 

444.8 

4.8 

2.7 

-142 

1 Apr 2011 (91) 

1507 

491.4 

5.4 

0.7 

-179 

27 Apr 2011 (117) 

1854 

475.5 

2.4 

19.0 

21 

14 May 2011 (134) 

2121 

549.0 

16.4 

15.7 

-85 

20 Jun 2011 (171) 

2349 

442.1 

4.8 

3.2 

47 

4 Jul 2011 (185) 

0630 

434.0 

2.5 

23.2 

-35 

9 Jul 2011 (190) 

0430 

477.7 

3.6 

22.0 

-65 

17 Jul 2011 (198) 

0120 

456.3 

3.3 

20.4 

-113 

25 Sep 2011 (268) 

1245 

516.5 

4.8 

0.3 

174 

1 Oct 2011 (274) 

0954 

485.7 

9.3 

22.8 

137 

6 Oct 2011 (279) 

0732 

500.2 

6.3 

21.7 

107 

12 Oct 2011 (285) 

0440 

527.9 

3.8 

20.5 

71 

25 Oct 2011 (298) 

1029 

481.2 

7.8 

18.3 

-10 


a From left to right the columns give the date and day of year (DOY), the 
time for which the Z p o~B p fit gives B p — 0, the fit intercept for Z p o, the 3cr 
standard error in Z p o, the MSO local time (LT) in hours, and the 
longitude of the B p = 0 point. 

internal fields contribute approximately equally to the total 
magnetic field, even at the planetary surface, the challenge in 
applying standard spherical harmonic analyses to derive Mer- 
cury’s internal field is substantial. Johnson et al. [2012] have 
discussed the derivation of an external model that is based on 
the parabolic magnetopause and magnetotail model of Alexeev 
et al. [2008, 2010] and for which the geometry and field 
intensity parameterizations are constrained by MESSENGER 
orbital observations. Here we present comparisons among 
candidate internal field models to which the external model 
derived by Johnson et al. [2012] has been added. 

3.1. Forward Model Residuals for Magnetically 
Quiet Cases 

[34] As Slavin et al. [2010] and numerous other authors 
have shown, magnetospheric dynamics at Mercury present 
strong variations in the observed magnetic field consistent 
with the presence of substantial scatter in the Z MSO versus B p 
data occasionally present (cf. Figures 3 and 8). We consider 
a subset of data selected to be as free of signatures of dynamics 
as possible, that is, for which the magnetic field data appear to 
vary smoothly throughout the magnetospheric transit and do 
not display prominent diamagnetic field depressions [e.g., 
Korth et al., 2011, 2012] or substorm signatures [Slavin et al., 
2010, 2012; Sundberg et al., 2012]. Moreover, we selected 
cases for which the standard error in the fit Z p0 was compar- 
atively small. In addition, events were chosen to span all 
planetary longitudes. The set of events selected for this com- 
parison is listed in Table 5. 

[35] To illustrate the constraints that the orbital data place 
on the internal magnetic field, we evaluated the predictions 
of five different internal field models along the trajectories 
of each magnetospheric transit corresponding to the events 
in Table 5. To each of these models we added the external 
field model of Alexeev et al. [2010] with the parameters 
determined by Johnson et al. [2012] from MESSENGER 
orbital data. The five internal field models are (1) the offset, 
axially aligned dipole of Johnson et al. [2012] with a 
southward directed moment of 190 nT-7?M and a northward 


offset of 0.196 R-m (he., that of Table 4); (2) an axially 
aligned, centered dipole with g- 10 = —190 nT; (3) an axially 
aligned, centered dipole with gi 0 = —300 nT, chosen to 
agree with the observed high-latitude field intensity; (4) an 
axially aligned, centered dipole-octupole field with g 10 = 
— 190 nT and g 30 = — 158 nT consistent with the g 30 inferred 
from the observed near-planet B p = 0 points; and (5) an 
axially aligned, centered dipole-octupole field with gi 0 = 

— 190 nT and g 30 = — 195 nT to illustrate the effects of an 
even stronger octupole contribution. 

[36] In Figure 13, we show a comparison of the B p and 
By observations for 26 March 2011 (DOY 85) from 
14:00:37 UTC (Z MSO = 1.0 R M ) to 14:23:45 UTC (Z MSO = 
—0.5 R m )- It is clear that the offset dipole represents the 
data well (Figure 13 (top)). The centered dipole with g 10 = 

— 1 90 nT yields a field that is too weak at higher Z MSO , and 
the centered dipole with g 10 = —300 nT yields a B p that is 
too broad and large at lower Z M so- The profiles of both 
centered dipole models are too broad, and only the offset 
dipole yields a field sufficiently focused at the higher Z mso . 
The B z component is also best fit by the offset dipole field. 

[37] Since higher-degree terms in a spherical harmonic 
expansion can also focus regions of increased field strength 
poleward, we consider the simplest such field, a dipole- 
octupole field (Figure 13 (bottom)), for two possible octupole 
contributions. Both dipole-octupole models yield a field that 
is too concentrated at the pole. Moreover, they both yield a 
much weaker field within ~0.4 R M of the equator than is 
observed. Both octupole models yield B p ~ 0 over a broad 
range in Z M so near the equator and a B z that is too weak near 
the equator. For both dipole-octupole models, dBJdZ ~ 0 
over a substantial region from Z MSO ~ 0 to 0.3 R M , which as 
noted above is not observed. Neither the B p nor B z behavior 
near the equator is consistent with the observations. The 
descrepancies are large, often reaching ~100 nT. 

[38] The same comparison between observations and the 
five internal field models is shown in Figure 14 for 25 October 
2011 (DOY 298) from 10:14:24 UTC (Z MSO = 1.0 R M ) to 
10:39: 13 UTC (Z MSO = —0.5 R M ). This event shows essentially 
the same results: only the offset dipole yields an appropriately 
focused field at higher Z MSO and a single clear magnetic 
equator. Only the residuals with the offset dipole are consis- 
tently low. All of the other models exhibit large systematic 
differences from the observations as functions of Z MSO , 

[39] The residuals at each Z M so were averaged over the 
orbits of Table 5, and the results are shown in Figure 15. 
Statistics given in Table 6 include the RMS residual, abso- 
lute residual, |<5|, and standard deviation of the residuals, cr-S, 
all averaged over Z MSO as indicated by 0. The centered 
dipole with gi 0 = —190 nT is not considered further since it 
is the most obviously wrong of the five models. The average 
standard deviations of the residuals, (cr-S), in B p , B z , and B T , 
are all in the range 4 to 1 1 nT. The offset dipole model yields 
(rms) values that are comparable to, or slightly larger than, 
the standard deviation in the residuals, indicating that any 
systematic signal exclusive of the offset dipole and external 
model is comparable to the (cr-S) values. 

[40] The average residuals versus Z MSO follow the same 
pattern as the examples. The centered dipole with g 10 = 
—300 nT gives B p = 50 to 100 nT more negative (positive 
residual) from Z M so = 0.2 to 0.8 R M , whereas B z is too large 
and positive (negative residual) near the equator, Z MSO < 
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Data 

g10=-300, RDZ=0 

g10=-190, RDZ=0.196 

— g10=-190, RDZ=0 




Figure 13. Comparison of observed magnetic field with different models for Mercury’s internal field for 
26 March 2011 (DOY 85) from 14:00:37 UTC (Z MS o = 1-0 R M ) to 14:23:45 UTC (Z MSO = -0.5 R M ). 
(left) Results for B p , and (right) results for B z . The observations are shown in black, (top) Results for three 
different dipole models: an offset dipole (red); a centered dipole that matches the high-latitude field (blue); 
and a centered dipole with the same moment as the offset dipole (green). RDZ is the dipole offset in R M . 
(bottom) Two centered dipole-octupole models, both with g w = — 190 nT: g 30 = — 158 nT consistent with the 
observed off-equator B p = 0 points (blue), and g° = — 1 95 nT (pink) to illustrate the effects of an even stronger 
octupole contribution. For reference: B p = 0 at 14:14:13 UTC when Z MSO = +0.180 R M (439 km) and 
B z = 65.8 nT; Z MSO = 0 at 14:16:44 UTC when B p = 24.1 nT and B z = 107.9 nT. The ascending-node 
MSO local time and planetary longitude were 2.7 h and — 142°E. 




Figure 14. Comparison of observed magnetic field with different models for Mercury’s internal field for 
25 October 2011 (DOY 298) from 10: 14:24 UTC (Z MS o = 1.0 R u ) to 10:39: 13 UTC (Z MS o = -0.5 R M ) in 
the same fomiat as Figure 13. For reference: B p = 0 at 10:28:36 UTC when Z MSO = +0.201 R M (490 km) 
and B z = 86.9 nT; Z M so = 0 at 10:31:37 UTC when B p = 21.1 nT and B z = 64.5 nT. The ascending-node 
MSO local time and planetary longitude were 18.3 h and — 10°E. 
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Figure 15. Average residuals of (top) B p and (bottom) B z 
components of the field for the orbits listed in Table 3 as func- 
tions of Z mso . The four models considered are the offset 
dipole model (red); the centered dipole model with a moment 
that matches the high northern latitude data (blue); the dipole- 
octupole model that gives an off-equatorial B p = 0 point where 
the near-planet crossings show the average B p = 0 point, g 30 = 
— 158 nT (solid black); and a dipole-octupole model with a 
larger octupole term, g 30 = —195 nT (dashed black). For the 
first two models, RDZ is the dipole offset in R M . The averages 
over Z MSO of the absolute residuals and of the standard devia- 
tions at each Z MSO are given in Table 6. In every case except 
the offset dipole model, the standard deviations in the residuals 
averaged over the cases in Table 5 are much smaller than the 
residuals. 

0.4 R m . Both of the dipole-octupole models underestimate 
B p for Z MSO = 0.2 to 0.8 R M by up to 70 nT, and then 
overestimate B p near Z MSO = 1 R M by 1 00 to 1 50 nT, whereas 
they underestimate B z for Z MSO <0.5 R M by ~50 nT and then 
overestimate B z near Z MSO = 0.8 to 0.9 R M by 60 to 100 nT. 
For these models, the (mis) and (|<5|) values are at least a 
factor of five greater than ( a-S ), indicating that the systematic 
residuals observed in the examples are similar for all of 
the events, that is, the discrepancies are similar at all 
longitudes. 

3.2. Upper Limits to Higher-Degree Structure 

[ 41 ] We use the residuals relative to the offset dipole field 
to assess any remaining signal that could be attributed to the 
axially aligned terms of degree 3 and 4 beyond that 


corresponding to the offset dipole (cf. Table 4). We denote 
these contributions by g 30 * and g 4 o*. Johnson et al. [2012] 
found that the residuals are best organized in MSO coordi- 
nates rather than body-fixed coordinates, indicating that most 
of the residual signal is due to external rather than internal 
fields. Thus, definitive estimation of nonzero higher degree 
and order terms must be deferred until the external model is 
refined to better than the level of residuals evident in 
Figure 15. Here we simply fit the residuals with a two- 
parameter axially symmetric field model, under the assump- 
tion that the signal can be attributed to a planetary origin. The 
parameters are the next two highest degrees in a spherical 
harmonic expansion of a zonal (axially symmetric) field: g 30 * 
and g 40 *. Because we attribute all of the residual signal to 
internal origin, and because we estimate only two parameters, 
the magnitudes that we find for g 30 * and g 40 * will almost 
certainly overestimate the amplitude of these terms in the 
internal field. Thus, the analysis here places an empirically 
derived upper bound on the magnitude of the additional 
contribution to the axial octupole term for qualitative com- 
parison with dynamo predictions. A similar approach has 
been recently used to analyze Saturn’s axisymmetric field 
\Cao et al, 201 1], 

[ 42 ] A grid search was perfonned to find the g 30 * and g 40 * 
coefficients for the best fit model to the residuals in B r of the 
offset dipole model relative to the observations for the orbits 
in Table 5. We restricted the data to portions of the orbit 
where the satellite altitude was less than 600 km, to ensure 
maximum sensitivity to the internal field. The resulting 
minimum RMS misfit was 7 nT for g 3 o*/gio = 0.04 and 
g 4 o*/gio = —0.07, from which we conclude that upper limits 
for the coefficients are |g 3 o*/gio|<0.04 and |g4o*/giol < 0-07, 
as given in Table 4. We emphasize that the g 40 * tern in 
particular likely contains power aliased from higher spheri- 
cal harmonic degrees and orders, and thus the results do not 
indicate the relative magnitude of g 30 * and g 40 *. The pri- 
mary result for our purposes is that |g 3 o*/g 2 ol is not larger 
than about 0.1, so that g 30 /g 2 o ~ 0.29 consistent with the 
offset dipole (cf. Tables 3 and 4). 

4. Summary and Discussion 

[ 43 ] The analyses of magnetic field equator signatures for 
equator crossings of the MESSENGER spacecraft on the 
descending node, at 1000 to 1500 km altitude, and ascending 
node, at 3500 to 5000 km altitude, give the same northward 
offset of the magnetic equator. The persistence of a north- 
ward offset in the distant altitude range confirms that the 
magnetic equator identification unambiguously determines 
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a Quantities in the table are defined as follows: <5 is the residual between the observation and the model for the specified internal 
field model and is a function of Z MSO ; () indicates averaging over Z MSO from 0 to 1 f? M ; (rms) is the root-mean-square misfit; 
(|<5J) is the average absolute value of 6; e-8 is the standard deviation in 8; and ( a-8 ) is its average over Z M so- 
The offset dipole model is that of Johnson et al. [2012]; moment = 190 nT -/?vi . Z-axis offset = +0.196 R M . 
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that the global field is best represented as a dipole offset to the 
north by 480 km. We note that the plasma sheet pressures are 
centered with respect to the magnetic equator [Korth et al., 
2012], providing additional confirmation of the offset. The 
most precise value for the offset is from the near-planet 
crossings, which give 479 ± 6 km, corresponding to g 2 J 
gio = 0.392 ± 0.010. This set of crossings yields a tilt of the 
offset dipole axis of not more than 0.8° from the planetary 
rotation axis. An axially aligned dipole-octupole field cannot 
explain the same offset of the near and far current sheet 
crossings from the planetary equator. Analysis of 13 mag- 
netically quiet orbits distributed approximately evenly in 
longitude reveals that the offset dipole representation with the 
external field of Johnson et al. [2012] yields residuals for 
these quiet orbits of ~10 nT, or ~2.5% of the maximum 
observed field in the northern polar regions. We find that g 30 
is dominated by the contribution from the offset dipole such 
thatg 30 /g 2 o ~ 0.3 andg 30 /gio ~ 0.1. This result is consistent 
with that obtained by Uno et al. [2009] from a regularized 
inversion of the Mariner 10 and MESSENGER flyby obser- 
vations. Our upper limits for the higher-degree axial coeffi- 
cients are also supported by inversions for high-latitude 
remanent crustal fields that yield g 30 * and g 40 * tenns that are 
less than 1% of the axial dipole term [Purucker et al., 2012], 

[44] The combination of high axial aligmnent, large value 
for |g 2 o/giol> an d low upper bound on ]g 3 o/giol presents a 
challenge to a number of dynamo models for Mercury’s 
magnetic field. As summarized in the introduction and 
reviewed by Stanley and Glatzmaier [2010] and Schubert 
and Soderlund [2011], many candidates for producing a 
weak field have been proposed. Among these models, those 
that predict power spectra dominated by odd hamionics, such 
as the induction feedback dynamo [Glassmeier et al., 2007a, 
2007b; Gomez-Perez and Solomon, 2010; Gomez-Perez and 
Wicht, 2010; Heyner et al., 2011a, 2011b] and the double 
iron-snow-layer model \ Vilim et al., 2010] are not compatible 
with either our large |g 2 o/giol or our low |g 3 o/g 2 ol ratios. 

[45] One mechanism that appears to account for the 
observed field invokes a non-convecting layer above a deep 
dynamo [Christensen, 2006; Christensen and Wicht, 2008], 
These models are variations on a suggestion originally pro- 
posed for Saturn [Stevenson, 1982]. Although the dynamo 
produces a field of high degree and order at the top of the 
convecting zone [cf. Heimpel et al., 2005], the field observed 
exterior to the core yields only lowest degree tenns because the 
higher-degree components of the field are sufficiently variable 
in time that they are highly attenuated as the field diffuses 
through the non-convecting, but electrically conductive, top 
layer of the core [Christensen, 2006; Christensen and Wicht, 
2008], In addition, differential rotation representing a thermal 
wind in a thick stable layer is important for preferentially 
attenuating the non-axisymmetric components of the low 
degree terms. However, the precise role of the stable layer in 
modifying the dynamo-produced field warrants further inves- 
tigation. For example, Stanley and Mohammadi [2008] found 
that thermal winds in a thin stable layer can destabilize the 
field morphology and lead to surface fields that need not be 
large in scale and can be either strong or weak. Moreover, 
relaxing the computational assumption of equal diffusivity 
for thermal and compositional buoyancy in the stable layer 
yields different surface field morphologies and strengths, 
depending on the core sulfur content [Manglik et al., 2010], 


For sulfur contents of 2%, the doubly diffusive convective 
regime for a thermally stably stratified layer produces fingers 
of compositionally light fluid that penetrate the stable layer, 
enhancing the poloidal field and resulting in a surface field that 
is an order of magnitude stronger than the observed field. 
Sulfur contents of 0.2% lead to weaker surface fields, but 
these fields are quadrupole-dominated rather than dipole- 
dominated. Thus these studies collectively indicate that 
understanding the dynamics of the stable layer is critical to 
predicting the observed surface field. 

[46] The recently determined gravity fields for Mercury 
[Smith et al., 2012] suggest that the solid outer shell of 
Mercury may be too dense to be composed solely of crustal 
and mantle silicates compatible with surface chemical 
observations and may contain an iron-rich lower layer, 
possibly in the form of solid Fe-S at the top of the core (S. A. 
Hauck II et al., The curious case of Mercury’s internal 
structure, submitted to Journal of Geophysical Research, 
2012). A solid layer, if sufficiently electrically conductive, 
might also attenuate the field structure at high degree and 
order. Possible core compositions consistent with the gravity 
results are also consistent with a stably stratified top layer in 
the liquid core such that the dynamo occurs in the con- 
vectively unstable lower portion of the fluid core [Smith 
et al., 2012; Hauck et al., submitted manuscript, 2012], 
This inference suggests that the model of Christensen and 
Wicht [2008], which follows the proposal for the axially 
aligned offset dipole field for Saturn [Stevenson, 1982], may 
warrant modification to allow for a solid layer overlying a 
stably stratified layer in the outermost core. 

[47] Measurement of the coefficients in the harmonic 
expansion of the internal magnetic field to higher degree and 
order may provide important clues to further discriminate 
among candidate dynamo scenarios, particularly in providing 
quantitative measures of higher-degree longitudinal struc- 
ture. The primary challenges to resolving structure at higher 
degree and order are limitations on current descriptions of the 
external field. Even at low altitudes the external field con- 
tributes as much to the observations as the internal field, so 
that errors of even a few percent in the external field will 
contribute residuals comparable to those found here. Johnson 
et al. [2012] determined a paraboloid model for the external 
field tuned to match orbital constraints for the magnetopause 
subsolar distance, tail current sheet distance and thickness, and 
magnetotail lobe flux. The residuals that remain after sub- 
tracting this external field and the offset dipole internal field 
are well organized in MSO coordinates but are not ordered in 
body-fixed coordinates. This finding implies that the primary 
remaining signal is from the external field, and Johnson et al. 
[2012] illustrated how a change in magnetopause distance 
can match some of the structure in the remaining residuals. 
R. M. Winslow et al. (Mercury’s magnetopause and bow 
shock from MESSENGER observations, submitted to Journal 
of Geophysical Research, 2012) examined the shape of 
Mercury’s magnetopause and bow shock and found that the 
parabolic magnetopause adopted in the external field model 
of Alexeev et al. [2008, 2010] is not an accurate description of 
the nightside magnetopause and that the magnetosphere is 
nearly cylindrical rather than parabolic on the nightside [cf. 
Shite et al., 1997], 

[48] Thus, it appears that the next major advance in char- 
acterizing the internal field depends on improving the 
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external field description and in particular the magnetopause 
shape. Nonetheless, the peculiarly simple global field implies 
that models yielding predominantly odd harmonics in the 
Gauss coefficients are at variance with the observations. 
Whereas the model that is most consistent qualitatively with 
the global field is a deep dynamo with an overlying non- 
convecting iron-rich layer [ Christensen , 2006; Christensen 
and Wicht, 2008], this model warrants further examination 
under less ideal assumptions, and continued efforts to match 
the observed structure with other mechanisms are needed to 
assess whether that structure is unique to the deep dynamo 
scenario. 
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